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Abstract We investigate the space-time evolution of the hydrodynamic particle-emitting sources with the fluctuating 


initial conditions generated by Heavy Ion Jet Interaction Generator (HIJING). In order to detect the event-by-event 


inhomogeneity of the sources, we examine the distribution of the error-inverse-weighted fluctuations, f, between the 


two-pion Bose-Einstein correlation functions of single and mixed events. We find that the distribution of f becomes 


wide for the fluctuating initial conditions. The large values of the distribution width and the root-mean-square fms are 


the signatures of the hydrodynamic particle-emitting sources with the fluctuating initial conditions. 
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1 Introduction 


Relativistic hydrodynamics has been widely used to 
describe the system evolution in high energy heavy ion 
collisions”, It provides the link between the initial 
and final states of the systems produced in the 
collisions. In general, the initial systems produced in 
relativistic heavy ion collisions are not spatially 
uniform, and there are event-by-event fluctuations of 
the initial quantities"). These initial fluctuations may 
affect the system evolution of space-time and lead to 
some changes of final particle observables relative to 
those associated with smoothed initial conditions *), 
Hanbury-Brown-Twiss (HBT) interferometry is 
a useful tool for probing the space-time structure of 
the particle-emitting sources in high energy heavy ion 


collisions” '*!, Previous studies indicate that the single 
event HBT correlation functions of the final identical 
pions exhibit event-by-event fluctuations in the 
[4 and Smoothed Particle 


Hydrodynamics (SPH) model!'*!, Investigating the 


granular source mode 


propagation of the initial fluctuations through system 
evolution and detecting their effects on the HBT 
measurements in relativistic heavy ion collisions are 
interesting issues. In this work, we use the Heavy Ion 


Jet Interaction Generator (HIJING)!"*! to generate the 
initial states of the systems event-by-event, at the 
energies of the Relativistic Heavy Ion Collider (RHIC) 
and the Large Hadron Collider (LHC). Then, we 
model the evolution of the particle-emitting sources 
with a (2+1) dimension relativistic hydrodynamics 
with the fluctuating initial conditions. Using the HBT 
technique based on event-by-event analysis''*!°), we 


investigate the effect of the fluctuating initial 


conditions on the space-time evolution of the 
found that the 


particle-emitting sources with HIJING fluctuating 


particle-emitting sources. It is 
initial conditions are inhomogeneous in space. The 
single-event two-pion correlation functions for the 
exhibit 
fluctuations. The large width of the distribution of the 
etror-inverse-weighted fluctuations, f, between the 


inhomogeneous sources event-by-event 


correlation functions of single and mixed events, and 
the large values of the root-mean-square fim; are the 
signatures of the particle-emitting sources with the 
fluctuating initial conditions. 


2 Hydrodynamic evolution of the systems 
with fluctuating initial conditions 


In the heavy ion collisions at RHIC and LHC energies, 
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the net baryon density of system is about zero because 
of the collision transparence. The ideal hydrodynamic 
description for the system with zero net baryon density 
is defined only by local energy and momentum 


[1,21 


conservation Under the l e of Bjorken 


longitudinal boost i invariance"! ”, the hydrodynamics in 
(3+1) dimension reduces to a (2+1) dimension 
hydrodynamics, and we need only to solve the 
(xy-plane) U8) 
Assuming that the system achieves local equilibrium 


transverse equations of motion 


at time tọ, we can use HIJING to construct the 


event-by-event initial energy density of the 


hydrodynamic evolving source at z=0 agl?70] 


Pa 1! 
&(t,,X, Y,Z=0 K) 
Foe i T h ano, 


[x= x E] Hy yE] 


20? 


(1) 


x exp 


where p,, is the transverse momentum of parton a, 


X,(T)) and y„(Tọ) are the transverse coordinates of 
the parton at t), o, is a transverse width parameter 
and K isa scale factor which can be adjusted to fit 
the experimental data of produced hadrons. 

Figure 1 shows the initial transverse 
distributions of two single events constructed with the 
HIJING at 7, = 0.6 fm/c for ./Syy =200 GeV Au+Au 
collisions with impact parameter b=4 fm and 8 fm, 
respectively. The parameters o, and K are taken 
to be 0.5fm and 0.8. One can see clearly that the 
transverse distribution of energy density for each event 
has large fluctuation, which leads to some "hot spots" 
in the initial sources. The number and maximum 
values of the spots decrease with the increase of 
impact parameter. 
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Fig.1 Initial transverse distributions of Pe density of two 
single events for Au+Au collisions at = 200 GeV, 
constructed with the HIJING at t) =0.6 fine, The panel (al) 
is for the impact parameter b=4 fm and the panel (a2) is for 
b=8 fm. The unit of energy density is GeV/fm’. 


To solve the hydrodynamic equations of 
motion, we also need the equation of state (EOS) 
which closes the hydrodynamic equations. In the 
calculations, we use the parameterized EOS named 
s95p-PCE”"! which combines the hadron resonance 
gas at low temperature and the lattice QCD results at 
high temperature. 


2.8 1263 
3 [ 


y/ fm 


-12 L111 1. -12 4 
“12-8-4 0 4 8 12 “12-8-4 0 4 8 12 


x / fm x / fm 


Fig.2 Transverse distributions of energy density of the two 
single events with b=4 fm (left panels) and 8 fm (right panels), 
at z=0 and t=3 fm/c [(b1) and (b2)], t=6 fm/c [(cl) and (c2)], 
and t=9 fim/c [(d1) and (d2)]. The initial conditions for panels 
(b1)-(d1) and (b2)-(d2) are the same as in Fig. 1(a1) and (a2), 
respectively. The unit of energy density is Gev/fm’. 

Once the initial conditions and EOS are 
determined, we can solve the hydrodynamic equations 
numerically. Figs.2(b1)-(d1) and (b2)-(d2) show the 
solutions of the transverse distributions of energy 
density at z=0 and t=3, 6, and 9 fm/c, for the two 
single-event sources with the HIJING initial 
conditions present in Figs.1(al) and (a2), respectively. 


[1,22] 


In the calculations, we ues the HLLE scheme and 


d"'273] to solve the 


the Sod's operator splitting metho 
hydrodynamic equations in the transverse plane z=0, 
and obtain the hydrodynamic solutions at z40 by the 
longitudinal boost invariance of  Bjorken 
hypothesis!'”'®!, The grid sizes for the HLLE is taken 
to be Ax=Ay=0.1 fm, and the time step is taken to 
be At=0.99-Ax for the Courant-Friedrichs-Lewy 
(CFL) criterion, At/Ax=v<1 [122] From Fig.2 it can 
be seen that the transverse distributions of the energy 
density have large event-by-event fluctuations in space. 
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These fluctuations are reserved even at the late stage 
of the evolution at t=9 fim/c. 


3 HBT signatures of 
particle-emitting sources 


inhomogeneous 


The HBT correlation function of identical pions is 
defined as the ratio of the two-pion momentum 
P(p,,p,) to the of the 
single-pion momentum distribution P(p,)P(p,) . 


distribution product 
Assuming that final identical pions are emitted at the 
space-time configuration characterized by a freeze-out 
temperature T,, we can generate the pion momentum 
P; (i=1,2) according to Bose-Einstein distribution, 
and then construct the single-event and mixed-event 


two-pion correlation functions“. 
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Fig.3 (Color online) Two-pion correlation functions for 
different single events (dashed lines) and mixed events (solid 
lines) with impact parameters b=4 fm (up panels) and b=8 fm 
(down panels). 


In Fig.3, we show our model-calculated 
two-pion correlation functions, C (qsides outs long)» for 
the single and mixed events with the HIJING initial 
Syy = 200 GeV AutAu collisions 
with b=4fm (up panels) and b=8 fm (down panels). 


conditions for 


Here the variables side, out, and Qiong are the 


39 


components of “side”, “out 
[24,25 


, and “long” of relative 
momentum of pion pair?*?°). The dash lines in each 
panel are the results for three different single events 
and the solid line is for the mixed event constructed 
with 80 single events. In the calculations, the 
freeze-out temperature is taken to be 120 MeV, and 
the total number of correlated pion-pairs in each event 
is Nyr=10°. One can see from Fig.3 that the correlation 
functions for the single events exhibit fluctuations 
relative to those for the mixed events, especially for 


the bigger impact parameter. This is because that the 


smaller number of the hot spots in the system with 
larger impact parameter may increase the source 
granularity”, and lead to larger fluctuations!®”!, 
However, in the usual mixed-event HBT analysis, 
these event-by-event fluctuations of single-event 
two-pion correlation functions are smoothed out. 

In order to observe the event-by-event 
fluctuations, we investigate the ratio of signal to noise 
of the fluctuations between the correlation functions of 


single and mixed events |C,(q,)—C,,(q;)| [14,15] 


|C, (qi) — Cn (Gi) | 
A| C,(q:)- Ca (4D I 


where i is the index of the component of relative 


f (qi) = (2) 


momentum q, and A|C,(q;)—C,,(q;)| is the error 
of the fluctuation. 
1000 
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Fig.4 (Color online) Distributions dN/df for 80 FIC and SIC 
events with impact parameters b = 4 fm (up panels) and b=8 
fm (down panels). N,,=10°. 

In Fig.4, we show the distributions of dN/df in 
the side, out, and long directions obtained from the 80 
events with impact parameter b=4 and 8, respectively. 
In calculations we take the width of the relative 
momentum q; bin to be 10 MeV/c and use the bins 
in the region 20<q, <200 MeV/c. The number of 
correlated pion pairs for each event is taken to be 
Nyn=10°. In Fig.4, the solid lines are the results for the 
fluctuating initial conditions (FIC). For comparison, 
the corresponding distributions for the events with the 
smoothed initial conditions (SIC), which are obtained 
by averaging over 400 random HIJING events, are 
shown with the dashed lines. One can see that the f 
distributions for FIC are wider than the corresponding 
results for SIC. The widths of the distributions for FIC 
increase with impact parameter. But the widths of the 
distributions for SIC are small and almost invariant 


with the increase of impact parameter. This is because 
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that the fluctuations of the single-event correlation 
functions for SIC are small always. 

For a limited Nan, we can reduce the number of 
analysis variable to decrease the noise in Eq.(2) and 
increase the ratio of signal to noise f, although it may 
lose some details. In Fig.5 we show the distributions 
of the f calculated with the variables of transverse 
relative momentum qp [panels (a) and (b)] and 
relative momentum q [panels (c) and (d)] of the pion 
pairs. The impact parameter b is 8 fm. One can see 
that the widths of the distributions of f for FIC increase 
with the number of the pion pairs in an event, Nw- 
Even for Non=10°, the widths for FIC are visibly larger 
than those for SIC. 


dN / df 


f (trans) f 


Fig.5 (Color online) Distributions dN/df for 80 FIC and SIC 
events with impact parameters b=8 fm , where f are 
calculated for variables qr [(a) and(b)] and q [(c) and (d)]. 

Next we examine quantitatively the widths of 
the distributions, dN/df, for the relativistic heavy ion 
collisions at RHIC and LHC energies. In Fig.6, we 
show the root-mean-square (RMS) of f as a function of 
Nyx calculated from the 40 simulated events for 
AutAu collisions at /syy = 200GeV and Pb+Pb 
collisions at A Shae = 2.76 TeV, respectively. 

From Fig.6 we can see that the values of fims 
for FIC increase with Nan. But for SIC the values of 
fims are almost independent of Ner. The reason is that 
for FIC the errors in Eq.(2) decrease with N, „ and 
for SIC both the differences |C,(q)—C,,(q)| and 
their errors decrease with Nawr. For the larger impact 
parameter, the fims for FIC increase more rapidly with 
Nx: because the number of the hot spots is smaller in 


this case. At LHC energy, we find that the hot-spot 
number in the source is larger than that at RHIC 
energy and with the same impact parameter. This leads 
to the smaller fims for FIC at LHC energy as compared 
to that at RHIC energy. 
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Fig.6 (Color online) Root-mean-square of f as a function of 
Non for the heavy ion collisions at RHIC and LHC energies. 

In experiments the number of correlated pion 
pairs in one event, N,,, is limited. For the central 
collisions at RHIC energy, the event multiplicity of 
identical pion, M, , is about several hundreds and N,, 
is about 10° (~M : / 2). However, at LHC energy, M 
is about several thousands and N,, may reach to 10°. 


T 


From Fig. 6 one can see that in these cases the results 
Of fims for FIC are visibly larger than those for SIC. 
of dN/df and fims for the 
inhomogeneous particle-emitting sources are hopefully 


The signatures 


to be detected in the heavy ion collisions at RHIC and 
LHC. 


4 Conclusion 


with HIJING 
conditions, we 


Using the model 


event-by-event fluctuating 


hydrodynamic 
initial 
investigate the space-time evolution of the 
particle-emitting sources in relativistic heavy ion 


collisions. The results indicate that the fluctuating 
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initial lead to 


inhomogeneous particle-emitting sources. For these 


conditions may event-by-event 
inhomogeneous sources the single-event two-pion 
exhibit 
However, in the usual HBT analyses performed for 


correlation functions large fluctuations. 
mixed events, these event-by-event fluctuations of the 
single-event correlation functions are smoothed out. 

In order to observe the fluctuations of the correlation 
functions, we investigate the distributions of f, the 
fluctuations between the two-pion correlation 
functions of single and mixed events with their 
error-inverse weights. We find that the widths of the 
distributions dN/df for FIC are much wider than those 
of the 


root-mean-square fms for FIC are larger. For FIC, fems 


for SIC. Correspondingly, the values 
increases with the impact parameter of collisions, 
because the smaller number of the hot spots in the 
system with larger impact parameter may lead to a 
larger source granularity. The values of fims for FIC 
increase with the number of correlated pion pairs, Nyx. 
However, the values of fims for SIC are almost 
independent of Nın. The large values of the distribution 
width and the root-mean-square fim; are the signatures 
of the particle-emitting source with the fluctuating 
initial conditions. Our calculations indicate that these 
signatures are hopefully to be detected in the heavy 
ion collisions at RHIC and LHC. 
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